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I. INTRODUCTION 


Nuclear matter in the spinodal region exhibits a universal behavior. When the uni¬ 
form matter enters into the spinodal region, small fluctuations grow rapidly and the matter 
undergoes a dynamical liquid-gas phase transformation. Nuclear multi-fragmentation pro¬ 
cesses are considered as a signature for such a dynamical phase transformation in nuclear 
matter [1-3]. Theoretical description of the spinodal dynamics requires an approach be¬ 
yond the mean-held theory by incorporating huctuation mechanisms [4-7]. The recently 
proposed stochastic mean-held approach (SMF) provides a useful framework for describing 
large amplitude nuclear collective motion by including one-body dissipation and huctuation 
mechanisms [8]. We carried out a number of investigations of the spinodal instabilities by 
employing the non-relativistic and relativistic framework of the SMF approach [9-13]. For 
a full description of the spinodal dynamics, we need to determine an ensemble of events by 
simulating the SMF equation with proper huctuations in the initial state. In our investiga¬ 
tion, instead of such simulations, we calculate the early growth of density huctuations in the 
linear response regime of the SMF approach. For this purpose, the equal time correlation 
fnnction of the density huctuations provides a very usefnl qnantity. For nnclear matter, 
in the linear response regime, it is possible to provide a nearly analytical treatment for 
the density correlation function. Indeed, the density correlation function exhibits the early 
stages of the dynamics of the liquid-gas phase transformation. In our earlier studies, we 
calculated the correlation fnnction in the pole approximation. As pointed ont by Bozek, the 
ehect of non-collective poles play an important role in the fnll description of the correlation 
function [14]. In a recent work, we calculated the density correlation function exactly by 
inclnding the ehects of non-collective poles in addition to the collective poles for charge sym¬ 
metric nuclear matter [15]. In this work, we present calculations of the density correlation 
fnnction for charge asymmetric nnclear matter in the semi-classical framework, i.e. in the 
long wavelength limit, of the SMF approach. In section 2, we briehy present description of 
the density correlation fnnctions including collective and non-collective modes in the linear 
response framework in the long wavelength limit. In section 3, we present calculations for 
the most nnstable initial density p = 0.3po at two different temperatures and two charge 
asymmetries of the matter. Conclusions are given in section 4. 
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II. DENSITY CORRELATION FUNCTIONS IN ASYMMETRIC NUCLEAR 
MATTER 


A detailed description of the correlation function of the density fluctuations for symmetric 
matter including the non-collective poles is presented in our recent publication [15]. Descrip¬ 
tion of the correlation function for the charge asymmetric system is given in a similar manner. 
Here, we give a brief description and refer to [9,15] for details. For charge asymmetric case, 
we need to consider the coupled SMF equations for neutrons and protons. The linearized 
equations for fluctuations of neutron and proton density matrices, Sp^(t) = p^if) — Pa are 
given by 


= [ha,6p^{t)] + [5U^{t),pa 


( 1 ) 


where pa is the density matrix of the initial state, ha denotes the mean-held Hamiltonian at 
the initial state and is the huctuating part of the mean-held for protons and neutrons 
a = p, n in the event labeled by A. Since, for inhnite matter the equilibrium state and the 
mean-held Hamiltonian are homogeneous, it is suitable to analyze these equations in plane 
wave representations. Then, the linear response equation becomes 

ih§-^<pi\5pl{t)\p2> ^2) 

= [^a(Pl) - £a(P2)] < Pl\5pl{t)\p2 > “ [pa{Pl) “ PaiVl)] < Pl\5U^{t)\p2 > 


The space Fourier transform of nucleon density huctuations 6p^{k,t) is related to the huc- 
tuations of the density matrix according to 

6pa{k,t) = f J <P + hk/2\6pa,s{t)\p- hk/2 >, (3) 

[My 

where the summation indicates sum over spin quantum number s =■]", i- We solve the linear 
response Eq. ([^ by employing the method of one-sided Fourier transformation in time [16- 
17], and obtain coupled algebraic equations for the Fourier transform of the local proton 
and neutron densities 6p^[k,u), 


( \ 


l + {k,0j) 

y [l + (fc, a;) 


^Pni^,^) + F7xn {k,u) Spp{k,u) 

6p^{k, u) + Py^Xp {k, uj) Sp^{k, uj) 


(4) 

S^ Ikoo) ) 


In deriving these coupled equations, we consider that the mean-held potential depends only 
on the local nucleon densities, = U (^Pp,p^^. The source term {k,oj^ is determined 
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by the matrix elements of fluctuations of the initial density matrix in spin-isospin 

channel as, 

d?p <f) + fik/2\5p^^g{Q)\'p-hk/2> 


S^Ak,oj)^Y. 


V ■ k — u 


(5) 


, . (27rh)3 

According to the basic postulate of the SMF approach, elements of the initial density matrix 
are uncorrelated Gaussian random numbers with zero mean values and with well-dehned 
variances. In semi-classical limit their variances are given by, 

< p + hk/2\6p^^^{0)\p — hfc/2 >< f — hP/2\6p^^^,{0)\f + Uk'/2 > 

= 5ab5ss'{!2.nh)^5{p- i^)6{hk - hP)fa{i^ (l - fa{p)) ■ (6) 


The factor Sab^ss' reflects the assumption that local density fluctuations in spin-isospin chan¬ 
nels are uncorrelated in the initial state. In Eq. (4), and Xa {k^uPj denote the zeroth 
order Landau parameters and the Lindhard functions for protons and neutrons, respectively. 
Landau parameters are dehned by the derivative of the mean-held potential energy with re¬ 
spect to the proton and neutron densities evaluated at the initial state, Fq^ = {dUb/dpa)^- 
The semi-classical expression of the Lindhard functions for neutrons and protons are given 
by 

(Fp p ■ hk/m dfa {p) 


Xaik,uj) = 2 


(7) 


{2ti hf hu - p-hk/m de 
Here the factor 2 comes from spin summation, and fa{p) is the Fermi-Dirac factor for protons 
or neutrons at the initial state. 

According to the method of one-sided Fourier transform [16-17], we can calculate time 
evolution of the local density huctuations by inverting Eq. Q and taking the inverse Fourier 
transform to give. 


Spaik,t) = -i 


r+o<^+i- dujG^{k,uj) 


iut 


( 8 ) 


/ — CX) + Z(T 27 r e{k,u) 

where the integral in the complex u is carried out along a contour that passes above all 
singularities of the integrand. The quantity G^{k,uj) in the integrand is given for neutrons 
and protons as 


{Pu) 

GilPoo) 


1 + {k, w)] (k, (u) - Fq^Xu {k, w) (k, w) 

Xn {k, w)] (fc, u) - F^Xp {k, ui) (k, u 


1 + K 


nn^ 

0 


(9) 


and e {k^ u?j denotes the susceptibility 


(fc, w) 


= 1 + i^n 


(k, u) + F^^xp {k, u) + [F^F^^ - F^^Fr] Xn {k, u) Xp {k, u) . (10) 
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As discussed in Ref. [15], it is possible to calculate the cu-integral in Eq. ([^, by choosing a 
contour as shown in Fig. 1 of Ref. [15]. There are collective poles determined by the roots 
of the dispersion relation, e{k,oj) = 0 —)■ a; = iiT^. The collective poles play important 
role in early growth of density fluctuations in the spinodal region. However, the collective 
poles alone do not provide a complete description of the growth of density fluctuations. The 
collective poles alone do not even satisfy the initial conditions as pointed in [14]. For a 
complete description, in addition to the collective poles, the effect of non-collective poles 
should be included into integral in Eq. ([^. By calculating the angular integration in the 
Lindhard functions, it is easy to see that the integrand has a logarithmic singularity. The 
integrand in complex cj-plane is multivalued, therefore entire real cn-axis is a branch cut. In 
order to calculate integral in Eq. ([^, we choose the contour C, as shown in Fig. 1 of Ref. 
[15]. We exclude the real cj-axis by drawing the contour from -|-cx) to the origin just above 
the real cn-axis, and after jumping from the hrst Riemann sheet to the second Riemann sheet 
at the origin, drawing contour just below the real (u-axis from origin to -|-cxd. Contour is 
completed with a large semi-circle and by jumping from the second Riemann sheet to the 
hrst one at origin. As a result, we hnd that the integral in Eq. ([^ can be expressed as 

= 6p^{P;k,t) + 5p^{C;k,t), (11) 


where the pole (P) contribution and cut (C) contribution are given by 

G^aik,±tTk) 


5Pa{P]k,t) = -Y, 

± de[k,oj)/doj\^=±ir^ 




and 


/•+“ doj 

Ga { k ,^ + i'n) 

G^{k,u-ip) 

1-00 271 

e{k^ oj + ip) 

e{k, OJ — ip) 


—iu)t 


( 12 ) 


(13) 


In order to investigate development of spinodal instabilities a very useful quantity is the 
equal time correlation functions of local density huctuations 

d?k 


= 5p^{f,t)5p^{r',t) = 


(2vr) 




(14) 


Here the local nucleon density huctuations Sp^{r', t)are determined by the Fourier transform 
of Sp^{k,t), and the spectral intensity aab{k,t)of the correlation function is dehned in terms 
of the variance of the Fourier transform of density huctuations according to 


aab{k, t){271 f 6{k - k') = 6p^{k,t)6p^{-kkt). 


(15) 


5 










In these expression, the bar indicates the average taken over the ensemble generated from 
the distribution of the initial fluctuations. 

We can calculate the spectral intensity aab{k, t) by evaluating the ensemble averages using 


Eqs. (12) and (13) for pole Sp^{P-, k,t) and cut 6p^{C; k,t) part of the Fourier transform 
density fluctuations and the Eq. @ for the initial fluctuations. As a result, the spectral 
intensity is expressed as. 


^ aab{PP',k,t) + 2aab[PC\ k,t) + aab(CC;k,t), 


(16) 


where the hrst and last term are due to pole and cut parts of the spectral intensity and the 
middle term denotes the mixed contribution. We briefly describe derivation of the analytical 
expressions of the various terms in the spectral intensity in the long wavelength limit in 
Appendix A. The total spectral intensity is obtained by summing over isospin components 


= <J„{k,() + 2t!^ik,t) + (T„n(k, t). 


( 17 ) 


The expression for the total correlation function of density fluctuations a{\r — which 

is summed over isospin components, is determined by using the total spectral density a{k, t) 
in Eq. (14). Using Eq. we can determine the initial condition of the spectral density 


a{k, 0), which must be equal to the right hand side of Eq. (18) at time f = 0. This leads to 
a non-trivial sum rule 


E 2 


p,n 


d?p 

(27rh)' 


-/a(^ (1 - /a(^) = o-pp{k, 0) + 2apn{k, 0) + ann{k, 0). 


(18) 


III. RESULTS OF CALCULATIONS 


In numerical calculations, we employ the effective Skyrme potential of Ref. [4] 

U^iPn, Pp) = A f + R f +c(^]Ta + l^^-DAp + D'Ap'Ta (19) 

\Po) \PoJ \poJ 2 dp Po 

where p = Pn + Pp and p’ = p^ — Pp are the total and relative densities, = +1 for neutrons 
and Ta = —1 for protons. The parameters A = —356.8 MeV, B = +303.9 MeV, a = 1/6 
and D = +130.0 MeV fm^ are adjusted to reproduce the saturation properties of nuclear 
matter: the binding energy eo = 15.7 MeV per nucleon, zero pressure at the saturation 
density po = 0.16 fm“^, compressibility coefficient K = 201 MeV and the surface energy 
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coefficient in the Weizsacker mass formula asurf = 18.6 MeV. Magnitude of the parameter 
D' = +34 MeV fm^ is close to the magnitude given by the SkM* [18]. The symmetry 
energy coefficient is C{p) = Ci — C' 2 (p/po)“ with Ci = +124.9 MeV and C 2 = +93.5 MeV. 
These values for the symmetry energy coefficient in the Weizsacker mass formula gives 
o-sym = £f{po)/^ + C'(Po)/2 = 36.9/3 + 31.4/2 = 28.0 MeV. We dehne the initial charge 
asymmetry as / = {pn — Pn) /P- 

Fig.l shows phase diagrams for different charge asymmetries corresponding to differ¬ 
ent wavelengths in temperature-density plane. These diagrams indicate the boundaries of 
spinodal unstable regions for different wavelengths, starting from upper most boundary for 
\ = oo. The critical temperatures corresponding to this effective Skyrme potential depends 
on the initial charge asymmetries I = 0.0, I = 0.4 and I = 0.8 are given by Tc = 15 MeV, 
Tc = 14 MeV and Tc = 10 MeV, respectively. The critical temperatures occur approxi¬ 
mately around the same initial density p = 0.3po- We notice that, consistent with earlier 
calculations [1], the spinodal instability region shrinks for increasing values of charge asym¬ 
metry. Neutron-rich matter with I = 0.8 and T = 1.0 MeV, approximately correspond to 
the structure of the crust of neutron stars. Under these conditions, the limiting spinodal 
boundary occurs at nucleon density around p = 0.55po, which is consistent with the result 
found in [19]. We pick a reference state located at the center of the spinodal region with 
an initial density p = 0.3po, and calculate the equal time correlation function of density 
fluctuations at two different temperatures T = 1.0 MeV and T = 5.0 MeV. In Figs. 2-3, we 
plot the total spectral intensity a{k, t) of correlation functions as a function of wave number 
k at time f = 40 fm/c for two different temperatures and three different initial charge asym¬ 
metries. At each initial charge asymmetry and temperatures, the upper limit of the wave 
number range /cmax is determined by the condition that the inverse growth rate of the mode 
vanishes, F^. = 0. Dashed, dash-dotted and solid lines indicate the result of calculations in 


Eq. (17) with pole contributions aab{PP', k, t) only, with cut contributions only and the total 
of all terms, respectively. The cut contributions include cut-cut contribution aabiCC] k,t) 
and the cross terms due to pole-cut parts 2aab{PC; k,t). From these hgures, similar to the 
charge symmetric matter, we make two important observations for each value of the charge 
asymmetry: Cut terms make an important negative contribution during the early phase of 
growth, hence slowing down the growth of instabilities. During later times, collective poles 
dominate the growth of density fluctuations, and the cut terms representing the effects of 
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non-collective poles do not grow in time, as discussed in earlier studies [1]. Second point 
that we note, both pole and cut contributions have divergent behavior with opposite sign, 
as wave numbers approach its upper limit, k —)• fcmax- These divergent behaviors cancel out 
each other to produce a nice regular behavior of the spectral intensity as a function of wave 
number. Fig. 4 shows the total spectral intensity of the correlation function a{k,t) for a 
typical charge asymmetry I = 0.4 as a function of wave number k at different times for 
initial temperatures T = 1.0 MeV and T = 5.0 MeV. We observe that the spectral intensity, 
which includes pole and cut contributions, clearly exhibits the growth of the initial density 
fluctuations. The initial value of the spectral intensity at time t = 0 indicated by the solid 
line is obtained by numerical calculations of the pole-pole terms, aab{PP', k,t = 0), cut-cut 
terms aab{CC] k,t = 0) and pole-cut terms 2aab{,PC] A;, t = 0). On the other hand, the solid 
circles on this line are calculated by the initial conditions located on the left hand side of 


Eq. (18). Nice agreement of both sides reflects the validity of highly non-trivial sum rule 


presented by Eq. (18). 

Fig. 5 shows the equal time correlation functions a{\f—r’\,t) of the total density fluctua¬ 
tions as a function distance of two space locations x = \r — r'\ at the initial density p = O.Spo 
and temperature T = 1.0 MeV for three different charge asymmetries. The density correla¬ 
tion functions are calculated with complete spectral intensities shown by solid lines in Fig. 
5, and the results are plotted at three different times t = 20, 30,40 fm/c. In Fig. 6, these cor¬ 
relation functions are plotted at temperature T = 5.0 MeV. The growth and organization of 
the correlation function describes the initial phase of the condensation mechanism. In order 
to recognize the condensation mechanism, it is useful to introduce the correlation length xc 
as the width of the correlation function at half maximum. The correlation length provides 
a measure for size of condensing droplets during growth of fluctuations. In the correlation 
volume AVc = dvrx^/S, the variance of local density fluctuations at time t is approximately 
given by a{xc,t). The number of nucleons in each correlation volume fluctuates with a 
dispersion AAc = AVc\l<j{xc-i t). As a result, the total nucleon number in each correlation 
volume fluctuates approximately within the range AAq — AAc < AA < AAq -|- AAc, where 
AAq = AVcpo denotes the number of nucleons at the initial uniform state. As seen from 
Figures 5 and 6, the correlation length is not very sensitive to the temperature and time evo¬ 
lution, but depends on the initial charge asymmetry. For example, at temperature T = 5.0 
MeV and the initial charge asymmetry I = 0.4, the correlation length is about xc = 3.0 





fm. In this case, the magnitude of the dispersion of density fluctuations at time t = 30 fm/c 
is about ^Jo^xc^ = 0.04 fm“^. Consequently, the number of nucleons in the correlation 
volume approximately fluctuates in the range of 1 < AA < 9. From this analysis, we observe 
that the spinodal decomposition indeed provides a dynamical mechanism for the liquid-gas 
phase transition. In the linear response approach, we can only recognize the early phase of 
the transition. In order to describe the full phase transition, we need to carry out long time 
simulations, which is not done in this study. 


IV. CONCLUSIONS 

As a continuation of a previous work, we investigate early growth of density fluctuations in 
charge asymmetric nuclear matter in the basis of the SMF approach. In the linear response 
framework, employing the method of one-sided Fourier transform, it is possible to carry 
out nearly analytical treatment of the correlation function of density fluctuations. The 
density correlation function provides very useful information about the initial phase of the 
liquid-gas phase transformation of the system. The method of one-sided Fourier transform 
involves a contour integration in the complex frequency plane. In earlier investigations, 
this contour integral is evaluated by keeping only effects of the collective poles associated 
with unstable collective modes. This approximation, although it describes the growth of 
density fluctuations in a satisfactory manner, has important drawbacks. It does not satisfy 
the initial conditions and furthermore leads to a divergent behavior as the wave numbers 
approach to their upper limits. In this work, we consider charge asymmetric nuclear matter 
in non-relativistic framework of the SMF approach. We calculate the correlation function of 
density fluctuations including the effects of collective poles and non-collective poles in terms 
of the cut contribution in the complex frequency plane. The cut contribution slows down the 
exponential growth of the pole contributions. It has also a divergent behavior with opposite 
sign, as wave numbers approach their upper limits. As a result divergent behaviors of pole 
and contributions cancel out each other to produce a nice regular behavior of the spectral 
intensity as a function of wave number. This allows us to have a complete description of 
the correlation function of density fluctuations in the linear response regime. Furthermore, 
exact calculations of the correlation functions satisfy the initial conditions expressed by a 
highly non-trivial sum rule. 
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Appendix A 


As seen in Eq. (17), there are three parts in the spectral intensity aab{k,t), the pole-pole 


part, the cut-cut part and the mixed cut-pole parts.The pole-pole part is 


aab{PP] k,t) = 


Kb 


de{k,u)/duj 


^ g-2rfct^ 


2E. 


ab 


ijj=iTi, 


de{k,uj)/duj 


(Al) 


- (jj=iTh 


where the quantities for neutrons and protons are given by 


Eln = [1 + FEK? + [Fo^Xn? /J, 


(A2) 


and 


El = [1 + Frxnf n + [ETXrf n, 


(A3) 


with 


= - (1 + E^^xv) FrXpI^ - (1 + F^Xn) F-^Xni;. 
~ d^p \{TkfT{p-k/my^ 


IT = 2 


(27r hy 


ifaip) (1 - faip)) ■ 


(A4) 


(AS) 


(rfc)2 -F (p- k/mj 

The pole-pole contributions have the same expressions as we had in our previous investigation 
in the semi-classical framework [16]. The cut-cut part has four terms. 


5„s(CC; k, t) = Alik, () + A+^(k, t) + Ajk, t) + t). 


(A6) 


The hrst and last terms come from the square of the hrst and second terms in Eq. (13), 


and involves double integrations over u and u'. The diagonal and off diagonal parts in the 
isospin space are 
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( ^nn{k,t) 


and 


duj r+°° du' 


1-00 271 J-oo 27r ca + o;' =F 2iT] 


/ITT 1 

1 ® i 

\ 1 

\ ^pn pp / 



(l)n.io^iri)+4>n {uj'^iri) 

e{uj^^iri)e{u!'^iri) 

tl)p{LLj^:iri)+(j>p{uj'^ip) 

e(u]^^ir})e{u}'^^ir]) 



Apnik,t) = - — 


+°° dw /■+°° dw' e 

— / - ( 

-oo 27T J-oo 271 U + u' ^ 2*7] V 


4>n. (cJ=F7Tj) + 0n 

e{ui^^ip)e{u]'^ip) 

(l>p{uj^ip)+(l>p{uj'^^ip) 

e{ui^^ip)e{u]'^ip) 


(A8) 


In these expressions and below, the symbol (8) denotes the matrix mnltiplication. The double 
integrals in A^fj{k,t) contain the principle value and the delta function contributions, which 
are identihed using the identity 1/ (ca + ca' =F = -P(l/a; + a;') ±77r(5(a; + a;'). The elements 
of W and V matrices are given by 


nn np 


pn pp 


and 



Vt \ / Fr.\f (i + Fry;) \ 

^ (1 + F„""xT) Fry;)' 

In these expressions and also below, we use the short hand notation yj 
x'a = Xa{k,oj' =F iv) the quantity 0a(i^ T iv) is dehned as 


(AlO) 
xyk,o^: ip), 


/ +00 


p-k/m — (uj^ irj) 


(All) 


There are contributions coming from the cross-terms in the square of Eq. (13), which 


are indicated by the second and third terms in Eq. (A6). These terms also involves double 
integrations over u and a;', and have similar structures as 


V ^ppik,t) 


r+°° du r+°° du' 

1—00 27r J—oo 27r oj -f- oj^ 


( \ 

I nn ^ np 

\ ^pn pp j 


^ <j>n ±ip) 

e{u)^iip)e{LLi' ±ip) 
4ip{u)Ofip)+<l>p {Lo'±ip) 
\ e{uiofip)e{iA}'±ip) 
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and 




dw dw' e 

—oo ‘2ti J—oo 27r u) -\- uj' 


pn nn 


(pn {L0±iri)+4in{uj'^iri) 
£{ui±iri)£{Lo' ^irf) 
pp {LO±iri)+pp 

£{(0±ir])£{Lo'^^iri) 


^ (A13) 


The integrand in in contrast to its appearance, is well-behaved func¬ 

tion, because when uj'=—u, the nominator is also zero therefore the ratio 
[0(a; =F iv) + irj)] / {oj + oj') becomes hnite. Here, the elements of W and V matri¬ 

ces are given by. 


nn * np 
* pn pp 


> + FZxp) (l + ITX 


(Frnpx' 


± 


{K’'rxpx'P 

(1 + F„””xJ) (l + F„””x' 



, (A14) 


and 


( ] 

\ ^ (l + FP-X'P ) ' 

In spectral density, the pole-cut contribution has also four terms 

dab{PC; k, t) = B+{k, t) + H+(fc, t) + B-^(k, t) + B-^(k, t). 
The isospin matrix elements of the first and the last term are given by. 


(A15) 


(A16) 


B]Pp(k, t) 




du 


de/duj\pj=xi^iT J-oo 27ra;=Ffr 


and 


/ X=F Y^] 

nn np 


Y^ X^ 

\ pn pp / 

\ 


(pn ip^'(^')~\'Pn 

I ('A17') 

PpiTiT')+4>p(‘^Tiv) * ’ ^ 
£{(^Tiv) 


B^ik,t) = 




r+°° du 


de/du\^=^ir J-oo 27ra;TirV 


X± Y± 


t'n k ')~\~4^n 

e{uj^ip) 

£(u^ip) 


(A18) 


Here , we use the short hand notation 0a(=Fir) = (f)a{k,u = T^T), ignore the label k in T^ 
for simplicity, and the elements of X and Y matrices are given by 


\ Y( AJ j \ (Frfxf^xp (l + •F’rxr'’) (1 + F^xD 
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and 


KT„ j ^ (l + Fo^xf) KVn j ' 


(A20) 


where Xa**" ~ Xa{k,oj = =F*r). The second and third terms involve an integral over u, and 
have similar structure as 


t) 


=F*e 


ipFi 


doj 6“*“^* 


and 


1 x^ ) 

nn np 


Y^ X^ 

\ pn pp / 

\ 


(/>n (=F»r)+(/>n (u]±ir]) 
£{LO±.iri) 

0p(=F»r)+(j!)p(LJifcit;) 

£{u)±.irj) 


(A21) 




p+oo g 


(9e/(9a;|a;=i|iir “'-oo 27r a; =f iT V 
The elements of X and Y matrices are given by 


XT, KT 


bn{TiY+4'n{‘^±iv) 

£{io±iri) 

i)p{^^ir)+<l>p{uj±iri) 

£{io±iri) 


(A22) 


X=F 

nn np 


Y^ 

pn pp 


and 


T Bq Xp ) "^0 Xp ) 

(Fon^rx^ 


{F7fx 

(l + F^^xf^) 


q=ir ± 

n An 

(1 + Frxi) 




(A23) 


[Y\^i FirxT^ (l + F^xt) t 
U>» / I (i + •f’rxj-'') Fs\i j ^ 


(A24) 
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FIG. 1: (Color online) Phase diagram in density-temperature plane for different wavelengths 


responding to the potential given by Eq. (19). 


cor- 
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FIG. 2: (Color online) Spectral intensity of the correlation function as a function of wave number 
at initial density p = O.Spo fm“^ at time t = 40 fm/c at temperature T = 1 MeV for three 
different charge asymmetries. Dotted, dashed-dotted and solid lines are results of pole, cut and 
total contributions, respectively. 
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FIG. 3: (Color online) Spectral intensity of the correlation function as a function of wave number 
at initial density p = 0.3po fm~^ at time t = 40 fm/c at temperature T = 5 MeV for three 
different charge asymmetries. Dotted, dashed-dotted and solid lines are results of pole, cut and 
total contributions, respectively. 
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FIG. 4: (Color online) Spectral intensity of the correlation function as a function of wave number 
at initial density p = O.Spo fm“^ and charge asymmetry I = 0.4 for different times at temperature 
T = 1 MeV (a) and T = 5 MeV (b). Dots on the solid lines at times t = 0 represent the initial 
conditions. 
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FIG. 5: (Color online)Density correlation function as a function of distance between two space loca¬ 
tion X = |r — f^l for initial density p = O.Spo fm”^ and different charge asymmetries at temperature 
T = 1 MeV at times t = 20, 30,40 fm/c. 
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FIG. 6: (Color online) Density correlation function as a function of distance between two space 
location x = \ for initial density p = 0.3/?o fm~^ and different charge asymmetries at tem¬ 

perature T = 5 MeV at times t = 20,30,40 fm/c. 
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